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1. Introduction 

Classical field theories that are invariant under scale and special conformal transformations gen¬ 
erally fail to retain these symmetries once quantized. Famously there is an anomaly, that is, the 
trace of the stress-energy tensor does not vanish, signaling a violation of invariance by rescalings. 
The exception consists of a class of quantum field theories for which the trace vanishes, known 
as conformal field theories (CFTs). When this happens not only is scale invariance restored, but 
the theory is also symmetric under the full group of conformal transformations [T]. This occurs 
at the fixed points of the renormalization group (RG) flow. 

It is also of interest to put the quantum field theory of interest on a curved background. When 
quantizing such a theory there are trace anomalies even at the fixed points of the RG flow of the 
corresponding flat-space theory [2] . These anomalies are given by a diffeomorphism-invariant local 
function involving derivatives of the metric. In d dimensions, there are a finite number of contri¬ 
butions of mass dimension d] for each there is a coefficient which is a function of the couplings. 
These coefficients are often of interest. Most notably, the coefficient a of the d-dimensional Euler 
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density is Cardy’s proposed extension [3] of the central charge c of two-dimensional CFTs, whose 
monotonicity properties under RG flow were understood by Zamolodchikov [3]. 

These coefficients, some of which are the central charges of the theory, are well-understood at 
the fixed points but are also defined along the RG flow. In two dimensions, a suitable extension 
of c away from fixed points may be defined, called c, which is a function of the couplings, so one 
may speak of their values along the RG flow in a sensible fashion. It is this quantity, c, that 
has the interesting properties that it decreases monotonically along RG flows and is stationary at 
fixed points where it takes the numerical value of the central charge c of the GFT corresponding 
to the fixed point. 

Given such remarkable properties of c, it is natural to ask whether such a quantity exists in 
the more physically interesting four-dimensional case. In fact, Weyl consistency conditions 
identify a quantity d in even spacetime dimensions that make it the one possible candidate for a 
generalization of Zamolodchikov’s c to higher dimensions. In four dimensions it was shown by Jack 
and Osborn [7] that this quantity is stationary at fixed points where it reduces to the coefficient 
a of the Euler term. Moreover, using pertnrbation theory they showed that this quantity is 
monotonically decreasing towards the IR. More specifically, they gave an equation for the RG 
flow of d that implies its monotonicity if a certain symmetric tensor, or “metric” in theory space 
parametrized by the couplings of the theory, is positive-definite. They then showed in an explicit 
perturbative calculation that this metric is in fact positive-definite for small conplings. More 
recently, positivity of this metric has been established in conformal perturbation theory BM- 

The extension of the quantity d to six dimensions was computed by a set of the current authors 
in [6], and furthermore was shown to have a natural definition in any even-dimensional spacetime 
as a consequence of the Weyl consistency conditions and the existence of a generalization of the 
Einstein tensor, along with a metric on the space of couplings that is analogous to that of Jack and 
Osborn and Zamolodchikov. This generalization of d is stationary at fixed points and rednces to 
a there. However, surprisingly, in m we showed by explicit computation in perturbation theory 
for a theory of scalars with a cubic self-coupling that the metric is negative-definite, and so d 
monotonically increases in the flow out of the trivial UV fixed point. Adding to this surprise, 
in [9] it was found that in a model with two-forms in six dimensions the metric is positive- 
definite. It seems that, even in perturbation theory, there is no straightforward generalization of 
the a-theorem in six dimensions, at least as envisioned in the cases so far. As explained in [9], this 
may be attributed to the fact that in six dimensions the trace anomaly on a conformal manifold 
defines three independent symmetric tensors on the space of couplings, only one of which satisfies 
positivity properties. This positive-definite tensor is, however, not the tensor that appears in the 
RG equation for a, and thus the monotonicity of its flow remains undetermined. Gontrary to 
this, in two and four dimensions there is a unique symmetric tensor with established positivity 
properties that also appears in the RG equation for c or d. 
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It is not known beyond perturbation theory whether flows of a in four dimensions are mono¬ 
tonic. However, there is another approach to the a-theorem that does not follow the previous lines 
of computation that uses unitarity of scattering processes in dilaton effective theories to establish 
positivity. Komargodski and Schwimmer have argued m without recourse to perturbation theory 
that the value of a on the UV fixed point is larger than that at the IR fixed point 0 A similar 
argument considered the same question in six dimensions m-, however it was not possible to 
reach a conclusion with the same methods as Komargodski and Schwimmer. Perhaps related in 
a general way to the difficulties encountered in m, we note that the (massless) scalar model 
with cubic interactions investigated in m has only a single Gaussian (trivial) fixed point within 
the domain of validity of the perturbative calculation, so the difference between the values of a 
in the UV and IR cannot be contemplated. Non-perturbative CFTs are known to exist in six 
dimensions, but since in addition to being non-perturbative they are non-Lagrangian CFTs, little 
is known about flows between them. 

As is clear from our discussion so far, the situation in six dimensions is significantly more 
complicated than that in two and four dimensions. We believe it would be useful to gain as much 
information as possible about the perturbative behavior of six-dimensional theories, beyond the 
computation of the quantity d. With this motivation in mind, we compute in this work, at two 
loops, the infinite part of the effective action and the trace anomaly in multiflavor cj)^ theory in 
six dimensions, including, for completeness of the analysis, the possibility that scale invariance is 
explicitly broken classically by a mass term. In addition to computing in a curved background 
(with a spacetime-dependent metric) we take the couplings to also have spacetime dependence. 
In effect this allows us to study the renormalization of correlators of operators that appear in 
the Lagrangian. With spacetime-dependent couplings counterterms proportional to derivatives of 
the couplings are required for finiteness. Correspondingly, the trace anomaly includes terms that 
contain derivatives not just of the metric but also of the couplings. The anomalies associated 
with these terms manifest themselves in the original model (with spacetime-independent metric 
and couplings) as coefficients of terms in the Greens’ functions of composite operators (including 
the stress-energy tensor and its trace). 

Given all these considerations, the focus of this paper is the Lagrangian 

of scalar fields (/>*, dehned on a six-dimensional manifold with metric 7 ^^^, where the repeated 
lowercase Latin flavor indices are to be summed over regardless of their position. This Lagrangian 
is of interest because it is the only general interacting theory that has classical scale invariance 


^This is a statement of the weak version of the a-theorem in four dimensions, where it is only established that 
auv > ctiR. The strong version of the a-theorem states that a decreases monotonically along the flow. The strongest 
version of the a-theorem states that this flow is a gradient flow |12| . 
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(for the appropriate choice of ^ij and zero rriij and hi) in six dimensionso Of course, one may 
consider theories that do not have Lagrangian descriptions in six dimensions, as mentioned above, 
but then the calculational methods of this paper are of no use and, typically, one must resort 
to holographic methods. With (|l.ip we may proceed in the old-fashioned ways of perturbation 
theory and reliably calculate the quantities of interest order by order in This is the starting 
point of this paper but first we must establish how such computations are performed on curved 
backgrounds. We should note that our results are reported here with the choice ^ij = ^6ij 
classically, with 6ij the Kronecker delta, as found from the general result ^ij = in d 

dimensions for conformal coupling of the scalar. 

The main computational method used in this work was developed and applied to various cases 
in four dimensions by Jack and Osborn in |15H18j . The main ingredients are the background field 
method and the heat kernel in dimensional regularization. In theory in six dimensions with 
a single scalar field and with spacetime-independent couplings, results for the two-loop effective 
action have been obtained in [IlHSI]; we have checked our results for some quantities against 
those listed in these references. 

The layout of this paper is as follows. In the next section we describe in detail the (perhaps 
unfamiliar but very powerful) computational method of Jack and Osborn. In section [3] we describe 
briefly the Weyl consistency conditions in order to make contact between our computations of the 
effective action and the a-theorem. In section H] we present our results for the infinite part of the 
effective action at two loops, and in section [5] we extract from those the two loop beta function 
and anomalous dimension. Finally, in section [ 6 ] we present results relevant to the a-theorem in 
six dimensions to three-loop order. Our conventions as well as various details and results needed 
for our computations are contained in three appendices. 

Note added in v2: The discussion in sections 3 and 6 contains errors that have been 
addressed in [ 22 ] . 


2. Method of calculation 

In this section we outline the method of calculation employed in this paper. For more details 
the reader is referred to HMS], where such computations have been thoroughly explained and 
demonstrated. Until section U] we assume for simplicity that no relevant parameters are present, 
for example rriij = 0 and hj = 0 in dni). 

In this work we will study quantum held theories dehned in spacetime dimension d = D — e, 
with D an integer, by a set of couplings g'^ and helds cj)'■ For our computations we will use 


^One may object that the theory is sick because of its potential, which is unbounded from below. However, 
within the context of perturbation theory, which is the scope of this paper, the ground state {(f>{x)) = 0 is stable to 
fluctuations of dfa;) |14| . 
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dimensional regularization and make explicit the mass dimension of the renormalized parameters 
via 




/ = 


( 2 . 1 ) 


for some numbers and 6 and where the index I labels the operators in the interaction La- 
grangian, i.e. I = (ijk) in (in])§ Though we start the perturbative calculations with g'^ and cj)', 
we will use (I2.1|) to express the resulting formulas in terms of the fields (j) and the dimensionless 
couplings . Then, with minimal subtraction, we have the bare parameters 


90 = + L\g )), 00 = , ( 2 . 2 ) 


with and — 1 containing just poles in e. In general is a matrix to account for the 
multiple number of fields in (ll.ip . The beta function and anomalous dimension are given by 

(3^ = fi——=—k^g^e + 13^ and 7 = de + —-— = de + 7 , (2.3) 

dg. dg 

respectively, where /3 and 7 are the quantum beta function and anomalous dimension respectively. 

Now, in quantum field theory in flat spacetime, wavefunction and coupling renormalization 
are enough to render finite correlation functions involving fundamental fields. When correlation 
functions involving composite operators are included, further counterterms are necessary. A conve¬ 
nient way to deal with these is by introducing sources for the composite operators, and including 
counterterms proportional to spacetime derivatives on those sources. For operators that appear in 
the Lagrangian it is enough to take their couplings as spacetime-dependent sources, g^ —>■ g\x), 
and introduce counterterms proportional to derivatives on g\x) [71123]. Finally, when a flat-space 
field theory is lifted to curved space with metric 7 ^,^, and the regularization procedure respects dif- 
feomorphism invariance, new divergences proportional to the curvatures defined from 7 ^jy appear, 
and thus further counterterms involving the curvatures are required for finiteness. 

In [7j a systematic treatment of such effects was undertaken, and a general expression for the 
Lagrangian in the presence of the sources 7 ^ 1 /(x) and g\x) was proposed, namely 


^0 = ^0 — g M + g , 


(2.4) 


where X-3% includes all field-independent counterterms, proportional only to curvatures and deriva¬ 
tives on g\x), and ^ = ^( 0 ) includes all field-dependent counterterms that also depend on 
curvatures and derivatives on g\x). jSfo is the bare Lagrangian, expressed in terms of g and 0 
with the use of () 2 . 2 I) . that contains terms that survive in flat space when the couplings are taken 
to be spacetime independent. It obeys the Callan-Symanzik equation 

= (2.5) 


®Note that the index carried by k of uni is not subject to the summation convention. 
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The RGE one finds from (j2.4p is 

(^'^7 - = ('S' ■ * + ('S^^T - ^ ■ P.6) 

which, by (j2.4p and the Callan-Symanzik equation (12.5p . requires 

(e-/3^^) A-.^ = /3 a-^, (2.7) 

and similarly for the terms, though there is an additional derivative with respect to the 

fields. As explained in [7] and we will review in the following, the terms /3 a • ^ defined by (12.711 
contribute, among others, to the trace anomaly of the theory in curved space. 

It is important to emphasize that in specific theories with possible relevant parameters like 
the RGE (12.61) is incomplete. For example, it does not correctly reproduce higher-order 
poles in higher-loop computations, even if the relevant parameters are set to zero in the classical 
Lagrangian. This issue has been analyzed in detail in [7] for four-dimensional theories, and also 
in [S] for (II.ip . While it does not affect our discussion below, it should be kept in mind. 


2.1. Background field method 

In this subsection we will give a brief overview of the background field method. We will present 
our expressions for the case of a single scalar field although the generalization to multiple fields 
and fields with spin is well-known. Our motivation for using the background field method is that 
it allows us to compute perturbatively counterterms like A • in (12.4p in a straightforward way. 

In the background field method one simply computes the effective action starting from 
which thus dictates the form of the counterterms. More specifically, we start by splitting the field 
(p into an arbitrary classical background part fib and a quantum fluctuation /, 


fi = fib + f ■ 


( 2 . 8 ) 


We can also introduce a source J, and obtain the effective action IE[</>b, J] (the generating func¬ 
tional of connected graphs with implicit and (x) dependence) after we integrate out 

/: 

^W[M = jDfe-So[4>]+fd‘^x^J{x)nx)^ j (-2.9) 

where 7 is the determinant of the metric which is not to be confused with the anomalous 
dimension 7 . 

To continue, let us denote by the action without any counterterms. Then, we expand 
in fluctuations. 





d’^xy/fi fMf + 5int[/], 


( 2 . 10 ) 


6 




where M = — with V the potential in Then, by expanding (j2.9p we find 

that, at the zeroth order, 

= ( 2 . 11 ) 

and at the one-loop order (a superscript in parentheses indicates the loop order), 

- ilndetM, (2.12) 

after we choose J appropriately in order to cancel terms linear in /, order by order in perturbation 
theory starting with ()2.10p . and subsequently perform in (j2.9p the Gaussian integral over /. Here, 
contains poles in e to cancel those in the — ^IndetM piece; in particular it contains the 
one-loop contributions to and L of (|2.2p . which are chosen to absorb the associated infinities 
coming from —^IndetM so that is finite. In addition, with the extension (12.41) it is clear 

from ()2.12p that also contains the one-loop contribution to that is given by the negative 
of the appropriate simple-pole part of — ^ In det M: 

I • i? C-(-ilndetM)P°'®. (2.13) 

Then, from (|2.7p and ()2.13l) we can evaluate Of course, — ^IndetM also contains field- 

dependent terms that require the counterterms for finiteness. 

At higher loops the interaction term S'int[/] in (I2.10p is considered and vacuum bubble dia¬ 
grams as well as diagrams with counterterm insertions are constructed^ The counterterms are 
of course fixed here by the previous loop order, i.e. by These diagrams can be evaluated in 
position space, using coincident limits of propagators according to the diagram topology. With 
these methods, which are explained thoroughly in the following, no loop integrations need to be 
performed. If we denote by the contribution of all such diagrams, we find 

IT(2) [,/,,] = -5f[</)fe] +=5^(2) _ (2.14) 

Again, hniteness of allows us to determine all counterterms in From the simple poles 

in A(2) . ^ it is again straightforward to evaluate ■ Si using the RGE (12.7p . Glearly these 
computations can be carried out order by order in perturbation theory. 

2.2. Heat kernel 

Using heat-kernel techniques the evaluation of (—^ Indet M)p°^® and higher loop poles may be 
accomplished. A pedagogical explanation of the method may be found in [25]; we will mainly 
follow the procedure as laid out in m, where it was used for single flavor theory without 


^As explained in [24], one of the advantages of computations done in the background field method is that only 
vacuum bubble diagrams need be considered order by order in perturbation theory. 
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spacetime-dependent couplings. A nice review of the heat kernel method and its applications can 
be found in [26] , 

The object of central importance in the heat kernel method is the propagator function in the 
presence of a background field as presented in section 12.11 It obeys the identity 

MikGkj{x,x') = 6ijS'^{x,x'), (2-15) 


where the indices are the flavor indices of the theory 5ij is the Kronecker delta, the d- 

dimensional biscalar delta function is defined by 

J d'^x'5^{x^x')4>{x') = (i){x), 7 '= 7 (x'), (2.16) 

and Mij is the elliptic differential operator, evaluated at the point x, alluded to in the previous 
section and defined by in our case of interest, having the general form 


Mij — -\- 


d‘^V{<p) 

d4d 


<l>=4‘b 


(2.17) 


The key to evaluating the determinant in the one-loop effective potential (I2.12p and the higher 
order diagrams, which involve integrals over products of Gij{x,x'), is to present Gij{x,x') in a 
way amenable to computation. 

The heat kernel provides such amenities. First, we define the heat kernel '^tj by the equation 


(^Sik^ +Mik^%j{x,x'-,t) = 0, (2.18) 

with 

'^ij{x,x'-jO) = 5ij6'^{x,x'). (2.19) 

Formally, by virtue of (I2.18p . the heat kernel may then be written 

^(t) = , ( 2 . 20 ) 

n 


with An the eigenvalues of Mij and the hats emphasizing that no particular basis of eigenstates 
I'i/’n) for the elliptic differential operator Mij need be chosen (however, the position basis will 
recover our calculations). Gij{x,x') may be then written as 

/•OO 

Gij{x,x') = / dt^ij{x,x';t). (2-21) 

Jo 

From the heat kernel ^ij{x,x']t) the one-loop effective action is obtained through the well- 
established zeta-function method (elaborated in, e.g., [T7] or m), which relates — ^ In det M to the 
heat kernel; to do so, an Ansatz for the form of the heat kernel must be given. This is suggested 








from the solution in flat space for the heat equation (|2.18p and was given by DeWitt [27] for a 
small t expansion: 




(A^+\d/2 ^ 2^an,tj[,X,X 


(47rt)'^/^ 


n=0 




( 2 . 22 ) 


with an,ij{x,x') the so-called Seeley-DeWitt coefficients and where a{x^x') is the biscalar distance- 
squared measure (called the geodetic interval by DeWitt), 

2 

, y(0) =x, y(l) = x', (2.23) 

0 * / 


a{x, x') = ^ 


dX 


'Yfj.v 


dy’' 
~aV 77V 


with y{X) a geodesic. AyM{x,x') is another biscalar, called the van Vleck-Morette determinant, 
that describes the spreading of geodesics from a point, defined by 

AvM(a:,x') = 7(x)"^/2 7(x')"^/^det • (2-24) 

We shall suppress the x, x' dependence of a, Avm, and an,ij henceforth. Now, the Ansatz (|2.22p 
obeys (12.1811 which yields the recursion relation 

nan,ij + dfj,a d^an,ij =with 9^ao,ij = 0 , (2.25) 


which allows us to compute the Seeley-DeWitt coefficients. 

With the asymptotic expansion of the propagator via the heat kernel established in (|2.22l) . its 
practicality in loop computations becomes evident. To elaborate on the comments above (12.2211 . at 
one loop one wishes to calculate the determinant in (|2.12p . This may accomplished by considering 
the so-called zeta function for the operator Mij, 

Cm{s) = J dtf~^jd‘^Xy^Wii{x,x-,t). (2.26) 

This function is useful to define because then the log of the determinant may be computed by 
differentiating it with respect to s and sending s to zero, which may be seen by considering the 
formal definition of in (| 2 . 20 ll : this yields 

— In det M = lim = f — [ d^xy^^ii{x,x]t). (2.27) 

^^0 ds Jo t J 

Given the formal definition in (I2.20jl , the value of In det M = In may computed with the 
equivalent of (I2.26P for with Tri^(t) = Explicitly evaluating (m as a function of 

s, differentiating with respect to s and then taking the limit as s —)• 0 reproduces the log of the 
determinant, formally, up to a minus sign. 

Actually, the equality in (12.2711 is only true up to the residue of a pole in s as s ^ 0 and 
eouation (12.271) is a bit misleadine at face value. Followine [T7]. limc^n C^(^) is of the form 

lim = lim ( [ dtt^~^ [ d'^x^'^ii{x,x]t) — — ) , (2.28) 

s— >-0 ds s— >-0 \Jq j s J 
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with P the residue of the integral inside the parentheses as s —0. However, in dimensional 
regularization this pole is displaced—this may be seen by noting the d dependence of the power 
series in t in (j2.22p . Hence, in dimensional regularization, P may be set to zero and we recover 
the e-dependent determinant 

dt C 

= — d'^Xy/^^ii{x,x-,t), (2.29) 

d=D-e Jo t J 

with D the integer dimension of spacetime, justifying the assertion in (I2.27p . Its pole, which is 
the main interest to us, may then be calculated with (|2.22l) and by noting the coincident limits 
therein, where x' x. \i D is even, then by expanding the series in (12.221) with d = D — e, it can 
be seen that the only piece that contains a pole in e as s ^ 0 is the (D/2)-th piece. Hence, the 
object of concern for the pole of the effective action is the coincident limit of the Seeley-DeWitt 
coefficient In six dimensions, in particular, we then have 

(-ilndetM);;°^" ^^g = ^7 / d^x^[as,ii]{x), (2.30) 

where the is inserted to preserve mass dimensions. The coincident limit of a 3 ^ii{x,x'), denoted 
in (j2.30l) by the brackets, may be found in appendix [Bl equation (IB.12p . and subsequently used 
to evaluate the one-loop counterterms of of (I2.12p . 

The task is then to extend the relatively graceful computation of the one-loop effective ac¬ 
tion, d la (j2.29D . to higher loop order. When using the heat kernel method in the context of 
the background field method, two-loop and higher-order contributions to the effective action are 
encompassed entirely within the calculation of vacuum bubble diagrams. These are then evalu¬ 
ated in coordinate space by integrations over the spacetime points involved in the loop diagram 
of the products of Green’s functions. It then becomes convenient, now specifying d = 6 — e, to 
express the Green’s function through the expansion (12.221) and (I2.2ip which, after performing the 
integration, yields 


(-Indet 


reg. 


1 - d(M 

= iim —;— 
ds 


Gij{x,x') = Go{x,x')ao^ij{x,x') + Gi{x,x') ai^ij{x,x') 

+ R2{x,x') a2,ij{x,x') -I- R^{x,x') a^^ij{x,x') + H{x,x '), 


(2.31) 


where the H(x,x') term does not contribute to UV divergences of the theory, i.e. do not have 
divergent behavior as x' ^ x§ The utility of this expansion is that it allows extraction of the 
poles of higher-loop diagrams almost by inspection, once the Gn{x,x') and Rn{x,x') are computed 
with (j2.21l) . For example, in the two-loop case, whose computation is detailed in section 01 the 


®It siiould be noted, however, that these still have IR-divergent behavior. Although it is not of interest in the 
calculation in this paper, it may be taken care of by considering the log of the ratio of the determinant of interest, 
equation (12.2911 . with determinant of the non-interacting operator = —SijV^. The ratio acts like a normalized 

version of (12.3011 and removes uninteresting IR divergences. See |15l section 1] for details. 
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coincident limit is necessary to evaluate the contribution there, so the computation boils down 
to the coincident limits of the Seeley-DeWitt coefficients, tabulated in the appendices, times 
coincident limits of the Gn{x,x') and which are easily computed with knowledge of 

the coincident limits of a{x,x') and tabulated in the appendices. Furthermore, 

although Gij{x,x') must be finite when a: 7 ^ x' in 6 or 6 — e dimensions, products of the Gn{x,x') 
and Rn{x,x') may have poles in e. In this two loop case the cubic product of Gij{x,x') is 
necessary, as evinced in Fig. [T] and (14.81) . and the various products of the pieces of (|2.3ip give 
rise to the poles computed in section [H 


2.3. Trace anomaly 


Now that the heat kernel method for the computation of the poles of the effective action has been 
established, we can proceed to the computation of the trace of the stress energy tensor defined by 
dni). To study the trace it is useful to promote the metric and couplings to spacetime dependent 


sources 




(2.32) 


and subsequently promote the action of the theory to be diffeomorphism invariant. Then we can 
define the quantum stress-energy tensor and finite composite operators by functionally differenti¬ 
ating with respect to these sources: 


T^u(x) = 2 - 




[Oiix)] = 


6S, 


0 


5g^{x) ’ 

where functional derivatives are defined in d spacetime dimensions by 

5\x, x '), J^ 9 '\x') = 6/S^{x, x'), 

with = ^{XjYj + XjYj). With these dehnitions it is easy to see that 

= eJ% + - (AcP) ■ ^So , 


(2.33) 


(2.34) 


(2.35) 


where A is the canonical scaling dimension of cj) and arises from variations of curvature- 
dependent terms in SYq. 

The trace anomaly may be viewed as the theory’s response to the local Weyl rescalings 


7^^(x) —)• (1 -|- 2fT(x)) 7^^(x), 5^(x) —)• {x) T- a{x)^\x). 


(2.36) 


The scalar ct(x) here is a variational parameter and should not be confused with the biscalar 
geodetic interval a{x,x') of the previous section. At the level of the generating functional we can 
implement these infinitesimal local Weyl transformations with the generators 


aW = 2 


j d'^Xy/^ 


ay 






4 = j d^^xViaP^j^. 


(2.37) 
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With these definitions 


A^W = -J , A^W = -j d'^x^a0^[Oj]). (2.38) 

Now, it is easy to see that the term e^o in (I2.35[l can be substituted with the use of (j2.6p . 
and so (|2.35l) can be written in the form 

- P^[Oi] D ^ - V^Z^), (2.39) 

where is the part of of (|2.35p that contains field-independent terms@ In (|2.39p we neglect 
field-dependent contributions besides those in j5^[Oi\. Equivalently, we can write 

J d^x^an-^Px-^ + J d'^x^ d^a (2.40) 

Terms in the right-hand side of (j2.39p have been computed in [7] for field theories in d = 4. In 
this work we will compute such terms for general multiflavor field theories in d = 6. As we 
just saw, these computations give results on the various terms that appear in the consistency 
conditions derived from ()2.40p [6]. 

Thus the relevant contributions to the trace of the stress energy tensor have their origin in 
the terms, which are, in turn, obtained from the heat kernel methods of the previous section. 
The Px ■ ^ terms are computed from the A • Si terms by (|2.7I) . and the Z^ terms are obtained 
from the Weyl variation - Si'). One can also change the basis so that terms in the variation 

dai^—X ■ Si') that appear in Z^ in one basis appear m fix - Si in another and vice-versa. 


3. Weyl consistency conditions 

The trace anomaly as presented in (I2.4n|] is useful because it allows very powerful statements about 
the structure of the theory along the renormalization group flow to be made. These statements 
arise from the Weyl consistency conditions, a specific example of the Wess-Zumino consistency 
conditions [28] that constrain the form of a quantum anomaly based upon the algebra of the 
anomalous symmetry group. 

Consider the two generators acting on the connected diagram generating functional W in (|2.40n . 
We may take the commutator of their actions on W for two different variational parameters a 
and a' and, because the Weyl rescalings are Abelian, obtain the Weyl consistency condition 

[AW-A^, AW-Af,]W = 0. (3.1) 

Now, the terms in IE, namely those coming from (|2.4p . have complicated transformations under 
(|2.37p and so (|3.ip imposes a set of non-trivial constraints and relations among these terms. 

®Note that the field-dependent part of is responsible for Oi —>■ [Oi], for the difference between d/dg^ and 5/Sg^ 
is a total derivative. 
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In particular, as argued in [5], must contain all terms that are diffeomorphism-invariant 

and, by simple power-counting, of mass dimension d that might arise in addition to the usual 
operators in in ( 12 .4|] from the promotion of the metric and couplings to spacetime-dependent 
sources as in ()2.32p . In two dimensions there are precisely three terms that fit the bill, with 
a single resulting consistency constraint. In four dimensions there are sixteen candidates, with 
seven independent consistency conditions. In six dimensions there are ninety five candidates with 
many independent consistency constraint equations. 

While the two- and four-dimensional cases are rather tractable and admit relatively simple 
interpretations of the consistency conditions, the six-dimensional case is significantly more complex. 
The f3\ ■ M and terms as well as the consistency conditions in six dimensions were categorized 
in [^. There, it was established that the most interesting consistency condition found in two and 
four dimensions, the consistency condition governing the flow of a certain o-function along the 
renormalization group flow, survives in six dimensions. The key point is to identify the analogs of 
the terms relevant to this consistency condition from two and four dimensions. There they involve 
the Euler density in the specific spacetime dimensionality and, in four dimensions, terms involving 
the Einstein tensor. In six dimensions, in fact, the coefficient a of the six-dimensional Euler term 
Eq is related to terms involving the generalization of the Einstein tensor, the so-called Lovelock 
tensor [29], in a way that is almosliJ completely analogous to the two- and four-dimensional 
settings. 

To be clear, in six dimensions (|2.40p takes the form | 6 ] 


65 » 30 „ 

(AW-Af)IT = / d^x^a{(3x-^)p + Yl / 

p=ld q=ld 

The terms of interest to the aforementioned consistency condition are contained therein 

(AW - A^)iy D J d^x^a(^-aEe - hL^ - 

+ [ d^Xy/jdi,a'H]dug^, 


(3.2) 


(3.3) 


where Li^ 3 , given in Appendix]^ are dimension-six curvature terms whose coefficients 61^3 vanish 
at fixed points, is the Lovelock curvature tensor, given in (IA.3p . and the coefficients T-L^j 
and are functions of the coupling constants. By varying each piece in (|3.3p with (|2.37p and 
applying the Weyl consistency conditions (|3.1I) one obtains the constraint equation 


did = + \{diV}j - djn])(5-^ , 


(3.4) 


^In six dimensions there are curvature terms with coefficients called “vanishing anomalies” whose only analog in 
lower dimensions is the lE term in four dimensions. At fixed points they vanish (hence the name) and the a in six 
dimensions is completely analogous to lower dimensions, however away from fixed points they modify a to 5 (to be 
discussed shortly) in a way distinct from two and four dimensions. 
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with 


(3.5) 


a = a + i 6 i - ^63 + . 

This equation is analogous to those found by Osborn in two and four dimensions [5]. 

By contracting with /3^ on each side of (|3.4I1 we arrive at six-dimensional equivalent of Zamolod- 
chikov’s theorem from two dimensions [ 1 ], 

= . (3.6) 

A similar relation was shown to hold in any even-dimensional spacetime in [ 6 ]. Thus, if we can 
compute Tijj in our theory and establish a definite sign, the monotonicity of the renormalization 
group flow of the theory can be established by way of (|3.6I1 . 

All that remains now is to determine the terms in the consistency conditions from the com¬ 
putation of the two-loop effective potential. The following section sets itself to this task. 


4. Poles of the effective action 

In this section we present our results for the pole part of the effective action up to two loops. 
For completeness we will also include here a mass term in our theory as well as a term linear in 
(j), i.e. we will take our Lagrangian to be given by 

+ {CijR + mij)4>i(pj) + hi4>i + ^gijk4>i4>j4>k , (4.1) 

where niij and hi have mass-dimension two and four respectively. Then, (12.4|) is modified to 

^Q=^Q-+ + (4.2) 

with ^ = ^{4>, g, 7 ) and ^ ^((/>, g^ m, 7 ). As in ( 12 .4p . in (|4.2p all quantities are bare quantities 
that may be written in terms of the renormalized quantities via (| 2 . 2 p . As stated in section \27I[ at 
one loop the effective action is related to the Seeley-DeWitt coefficient 0,012,ij with D the even 
integer spacetime dimension. To wit, in the theory of (14.11) in six dimensions, as in (I2.30p . 

(-4lndetM)P°i'^ = J d'^x^ [a,,ii]ix). (4.3) 

This result produces all terms in (14.211 . Thus, at one loop, using the result from appendix iBl 
discarding total derivatives, and choosing ^ij = (4 — 6ij from here on|®| we can isolate each 


®The e-dependent portion is necessary to maintain, in the absence of iriij and hi, classical conformal invariance of 
dm in 6 — e dimensions. Although inconsequential at one loop, this e-dependence is crucial at two loops. This has 
also been seen in higher-loop computations in four dimensions in [30]. We thank Hugh Osborn for bringing this issue 
to our attention. 
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piece of the one-loop effective action according to (I4.2p and (|2.2I) . We finco 

^ ^ (“9072-^6 + 540^1 “ 3(m'^2 “ > (4.4) 

where is the number of scalar fields (p. The field-dependent counterterms at one loop are given 
by 

~ 647 r^ 6 9ikl9jkl i ^ijk ~ 7 647 r^ {dimndjmpdknp ~ i2^9ijl9kmn9lmn + permutations)) , 

(4.5) 

and 

+ Jo R9iki9jki(pi(pj + \d^gikid^9jki(pi(pj + gikid^gjkidfj^pipj) . (4.6) 

Finally, the mass-dependent counterterms are 
^ ~ ~7 647r^ 2 (M ^~b 

T T5(('^ ~ 9ijk4’k T gikl9jkm4^l4^m T 'kk^ij'^klik 9jkl4^l (4.7) 

-h Id^mijd^niij + \mijmjkmki) ■ 

At two loops in the background field method we must compute the relevant vacuum bubble 
diagrams. Thus we are led to consider the diagrams in Fig. [H where in the diagram on the 
right 0 denotes the one-loop counterterm. Note that these are graphs in position space, and 

X x' 

Fig. 1 : The diagrams that need to be considered at the two-loop level. 

that short distance singularities arise here from the coincident limit of products of position-space 
propagators. In particular, the left graph in Fig. [T] is given by 

Vy 9ijk{x)9imn{x')Gii{x,x')Gjm{x,x')Gknix,x') , (4.8) 

with the factor of a symmetry factor from interchanging the Green’s functions. The evaluation 
of this integral is straightforward, once the divergent parts of the products of the Gn{x,x') from 
(j2.31l) are known; as these are listed in appendix [Cl we will not explicitly show the intermediate 
details of the cube of the propagator in ([4.80 and simply include their contributions to the 
counterterms in (j4.2|i in the results listed below. 

The graph on the right in Fig. [1] is slightly more complicated than the expression listed 
in m because of the spacetime dependence of the couplings. The counterterm insertion to the 

®This agrees with the results of |31| . where similar computations were performed for fermions and two-forms. 
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propagator in that diagram is simply obtained from our one-loop results. More specifically, using 
our results (14.6p and (021), or, equivalently, the coincident limit of a 3 ^ii{x,x') in (IB.12p . we obtain 



1 1 




d X d X 's/^ 5 (x, X ) 2 (nZmn^n 4“ 4“ (C/m g 


4“ 12^^''^ (^9ikl(.^)9jkm{^ ')Gij(^X,X . 

(4.9) 

Here we have specified the spacetime point at which any derivatives are to be taken. 

The UV divergences of these graphs are obtained from the divergences that arise from the 
products of propagators, discussed at the end of section 12.21 We use the results listed in ap¬ 
pendix [U1 along with prudent integrations by parts, to reduce both (14.81) and (14.9p to poles in e 
and coincident limits of the Seeley-DeWitt coefficients, which ultimately give explicit expression 
for the and pieces of (14.2p . In particular, for the counterterm 

graph, we finJ^ 

1 /i“ 



d ^ 2 9ikl9jkm (^9lmn4^n 4“ IXLlm 4“ (C/m 3g<J/m) 


(4.10) 


(647r3)2 

12 9ikl ■ 2 [c12,y] ( 3 /) 

4“ 12 9iki9jki (^[^ ( 3 /) (6 e) [ 03 ^ jj] ( 3 /)) 

Evaluting these graphs with the coincident limits listed in appendices [B] and [Cl we generate 
the terms as listed in ()4.2p at the two loop level. For the curvature and derivatives on couplings 
terms, at two loops discarding total derivatives we have 

4- 9ijk dfiX7^9ijk ~ ^dygijk^'^gijk 4- R^^g-ijk'^^9ijk) 

- + 30 - 10 - 3H^^)d^g,jkdugijk 

+ {E, + §-,F+^^R^- i^V^R) d^g,,kd,g^jk) ■ 

This result has been given in a conformally-covariant basis in [9]. Field-dependent counterterms 
at two loops are given by 

^ij ~ ~7 ( 647 r^)^ TS {,9ikl9jmn9kmp9lnp ^ 9ikl9jkm9lnp9mnp} , (^-1^) 

^ijk ~ ~7 ( 647 r^)-^ (^'g{9ilm9jnpgkqr9lnq9mpr 4“ 2 9ilm9jln9kpq9mprgnqr J 2 9ilm9jln9kmp9nqrgpqr) 

^i9ijl9kmn9lpq9mpr9nqr ^9ijl9kmn9lmp9nqr9pqr 4“ permutations)) , 

(4.13) 


(4.11) 


^'^This equation may be compared with (3.14b) of |19| . Note that our result reduces to this equation in the single 
spacetime-independent coupling case. 
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which are relevant for the two-loop anomalous dimension and beta function, and 

^ ^ ~ “7 ( 647 r^)^ 18 (( 288 O ^9ijk9klm9jlm 43 9ijkd^gjmndu9kmn 

^ ^9ijk9jlm^ 9klm 2400 Sijk 9^9jlm d^9klm 

9jlm^ 9klm 9ijk^^d 9jlm^9klm 

'^9ijkd^^ 9jlmdfi9klm iq 9ijk9jlm^ ^ gklrnj^ki 

i.'^2j(^^9ikl9jmn9kmp9lnp I^9ikl9jkm9lnp9mnp 

29klm9knpd^ 9ilndp9jmp 9klm9klnd^ 9impdp9jnp 
2 9ikl9kmnd^9jlp9p9mnp '^9ikl9km,nd^9jmpdpL9lnp 
^ 9ikl9mnpd^9jkmdp9lnp 'g9ikl9jmnd^9kmpdfj,9lnp 

IQ 9ikl9jkmd^9lnpd^9mnp^'Pi4'j 

“1“ (^9jkl9kmn9l mp d^9inp + ^9ikl9j mnQkmp 9^9lnp 

9jkm9mnpd^ (^giklQlnp) 2 9ikl9jkrn.9lnp9^ 9mnp) 4^i9^(pj^ . 

The mass-dependent counterterms are given by 

^ 7 (647r^)^ TS ((2880 ^9ikl9jkl 48 d^gikldy9jkl 

“1“ '^^9ikl^ 9jkl T 2400 ^9^ 9ikl 9^9jkl 
- ^^‘^9ikN‘^9jki - ^^^9''gikiV^9ygjki - \9^giki9^V‘^gjki 
T 10'^^ 9iki9^gjki9i/ Rgiki9^gjki9^ T yq 9iki9^^ gjki9^ 

T IQQ R9ikl9jkm9lmn4^n R9ikl9jmn9kmn4^l 

4 9ikl^ 9jkm9lmn4^n T '^9ikl^ 9jmn9kmn^l 

T 2 9ikl9j mn ^ 9kmn 4^1 

^9^9ikl9^gjkm9lmn4^n T q 9ikl9^gjmn9^gmnk4^l 

‘^9ikl9^ gjkm9lmn4^n9^ + 9ikl9^ gjmn9kmn4^l9^ 

T 2 9ikl9j mn d^9k mn 4>idi, 

~ 9ikl9jkm9lmn(pn^ “1“ gikWjmngkmn't’l^ 

4 9ikl9j mn9kmp9lnq4^p4^q T g9ikl9jkn9lpq9npr4^q4^r 
+ Y,9ikl9j mn9kmp9npq4^l4^q g9ikl9jmn 9kpq9mnp 4^l4^q 

“t“ ( 200 0^ ^9imn9kmn^jl 

T 4,9^gikm9^gjlm Y^9^gimn9^gkmn^jl 


(4.14) 
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“ 1 “ '^9imn9jkp9lmp4’n 'g9ikm9jln9mnp'Pp 9ikm9jnp9lnp4^m 
“ 1 “ ^9imn9kmp9npq4^q^jl g 9imn9kpq9mnp'i^q^jl)'^ij'^kl 
(^4 9ikmd^9jlm 2 9imnd^9kmn^jl '^9^9imn9kmn^jl^'^ijdp^kl 

“ 1 “ (^2 9ikm9jlm 9imn9kmn^jl) d^^ij dfi^^kl 

(^4 9ikm9jln 9ikp9lmp^jn fg 9ipq9kpq^jm^ln)'^ij'IT^kl'IT^mnj • 


Although these terms are unsightly, they allow us to calculate the quantities of interest—they 
give us the complete, general trace anomaly on a curved background with spacetime dependent 
marginal sources { 9 ijk{x) and 7 ^jy(x)) via equation (|2.39p . Each of the terms presented in this 
section yields the relevant beta functions in the first part of (I2.39P via equation (12.71) . The second 
set of terms in ()2.39p . called are obtained from a Weyl variation of the \-M terms, as seen 

in (j2.40p . Since the Weyl variation is non-trivial, we report the one- and two-loop contribution 
to here, since it is required for the identification of terms necessary to the computation of a 
in section El 

At the one-loop level there are no contributions from the Weyl variation of (|4.4p and so 
= 0. At two loops the Weyl variation of (14.111) yields 

^ ~ (647ri>)^ 1620 (^^4,9ijk 9^ 9ijk T 4 ^ ^ 9ijkd^ 9ijk T 400 '^ 9ijkd^ 9ijk R9ijkd^ 9ijk 

8-^1 9ijkdp9ijk ~4 ^2 9ijkdu9ijk T -^3 9ijk^u9ijk 

T ' 4^4 9ijkdi/9ijk T 8^5 9ijkdi/9ijk 
+ lV>^R’'P9ijk^udpgijk + %V^R''Pd,gijkdpgijk - R^^^,g^jk^pg^jk 

+ ^d'^Rd^gijkd^gijk + ^d^^Rdt,gijkd''gijk 
+ lR''^9ijk'^^^vdpgijk + ^R^"'gijkdu^'^gijk - ^R^'^d^^gijk'^udpgijk 

- iR^'^dugijk^^gijk + ^R^Pd^'gijkVudpgijk + R!'^d,gijkV^dpgijk 

- ^Rgijkd^^'^gijk + ^Rd^gijk'^'^gijk - ^R'^^dygi^kd'"g^u 

+ iR'^^'^'^dpgijkVyd^gijk 

- vdpgijkV^S/'"dPgijk -h gijkdyS/'^gijk + ^^'^gijkd^S/^gijk 

- ^dygijkV^^d’'V‘^gijk + 

+ ^9ijkd>"'^‘^V^9ijk'j ■ 

(4.16) 

It should be noted that the basis reported here is not identical to the terms reported as in 
p3.2p . which refers to the basis used in [ 6 ] written by some of the authors of the present work. 
However, as that basis is complete, the terms in (j4.16p may be written in the basis with 
repeated and judicious integrations by parts. For the purposes of the calculations in section [ 6 ] of 
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this paper, we did not find the basis referred to in ()3.2I) useful and were able to identify those 
terms required in equations (I6.6al) and (|6.6bh from the current presentation in (|4.16l) . 

From equations (I4.4jl to (14.7p and (14.lip to (14.1511 we can extract the beta functions for the 
couplings and masses, the anomalous dimensions of the fields, and, perhaps most importantly, the 
quantity a (and d), as described in section [31 which is the analog of Zamolodchikov’s celebrated 
c. 


5. Beta functions and anomalous dimensions 


As we have seen using background field and heat kernel methods the computation of Zij and Lijk 
is easily done in position space and does not require the calculation of any integrals. With our 
results (14.511 . ()4.12ll and (|4.13ll we can now compute the anomalous dimension jij of (pi and the 
beta function f3ijk to two-loop order. 

The anomalous dimension is defined by 


7 




dt 


t = - ln{fi/fio), 


where the RG time t is defined to increase as we flow to the IR. At one loop we find 

1 1 


7(1) = 


647r^ 12 

where we use the diagram to denote the corresponding contraction of the couplings, i.e. 


diklQjkl 


The two-loop anomalous dimension is 


^( 2 ) = 


1 


1 


11 

2i 


(647r3)2 18 

For the case of a single field cp our results (j5.2p and (j5.4p reduce to the results of 
also [I9H2I1133]). 

The beta function is defined by 


Of i dg dg 


At one loop we find 


i 

647r3 y 



1 

12 



(5.1) 

(5.2) 

(5.3) 

(5.4) 
(see 

(5.5) 

(5.6) 


where permutations of the free indices in the wavefunction-renormalization correction are under¬ 
stood, i.e. 

= gijigimngkmn + permutations . (5.7) 
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Eq. (15.611 reproduces the result of [32] (see also [T9H2n[33] l in the case of a single field cj)- In that 
case has a negative sign, and hence the corresponding theory is asymptotically-free. The 
two-loop beta function is 



(5.8) 

The first contribution to (15.8p is non-planar. For the seemingly asymmetric vertex corrections in 
(|5.8p (second and third term) a symmetrization is understood; for example, 



represents 



(5.9) 


In the single-field case (|5.8p reproduces the result of [32] (see also [T^2lj^^h . which, just like 
is also negative. 

The results presented here for the anomalous dimension and the beta function to two loops 
are found to agree with the results of [M], and can also be fully extracted from [35j . 


6. The metric in coupling space and the a-anomaly 

In section [3l and in particular in equation (13.61) , it was made apparent that there is an important 
piece of the /3\-M terms, called Tijj in this paper and in |6|@ that manifests itself as the coefficient 
of contact terms of certain correlation functions of the operators of the theory in flat spacetime 
and spacetime-independent gijk{x) and mij{x). This metric is important because it controls the 
behavior of a (or really, o) along the renormalization group flow; given the outstanding importance 
given to a (or its analogs) in two and four dimensions for its central role in characterizing quantum 
field theories there, its behavior in six dimensions gives insight into the universal features of 
quantum field theories in any dimension, possibly beyond the conventional Lagrangian description 
so ubiquitous in our understanding today. 

In [TU] a perturbative computation of the theory defined in the Lagrangian formalism by (II. ip 
yielded a surprising result for the value of 'H]j at the two loop level. One of the main purposes 
of this work is to give the details of that computation, along with other interesting results from 
the computation of the effective action. 

In order to compute we must first identify the corresponding piece in so that we may 
use (12. 7p to obtain the quantity of interest. There is no candidate in the one loop computation, 
Thus, we must look for a two loop contribution in where we do indeed find a 


There is a typo in the relevant equation in |19 |. 

was called xu in the two- and four-dimensional cases of [5] and in the six-dimensional case of |10| . 
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candidate. We see from (|4.1ip that the relevant piece is 


A( 2 ) — 


1 1 


1 


dfiQijk du9ijk 


( 6 . 1 ) 


e (647r3)2 1296 O' 

From (j6.ip and (I2.40p we can immediately match to the term \T-L\jd^g^dvg'^H^'^ in and 

extract 


^ 1 ( 2 ) _ 
^ij — ' 


1 


1 


:^IJ 


(647r3)2 3240"^'^ ’ 

where, as section [2l we use notation of [6] and denote I = (ijk). Furthermore, performing a Weyl 
variation of dSI]) we find 

1 1 


-M) D - 


e (647r3)2 5430 


dugijk dfj,a, 


and so 


1 


1 


^ 1 ( 2 ) ^ _ 

^ (647r3)2 12960 


91 : 


(6.3) 


(6.4) 


as also seen in (I4.16p . The result (j6.2p is unambiguous and scheme-independent. As we observe, 
the leading, two-loop contribution to the metric is negative, and so the consistency condition (13.4p 
and its consequence (13.6p cannot possibly lead to a strong a-theorem for a. 

Now, our theory has only the Gaussian fixed point in perturbation theory. Non-perturbatively 
there may be a non-trivial fixed point, but our results (16.2p and p6.4p cannot be used beyond per¬ 
turbation theory. Nevertheless, as long as the flow of our theory can be described perturbatively, 
the quantity a is monotonically increasing. 

Another use of the consistency conditions is the evaluation of some quantities at higher loop 
orders. Regarding a, for example, we can use (13.4p with the results (j5.6h . (j6.2p . and (16.4h to 
obtain the three-loop contribution to d, 


a(3) = 


(647r3)3 77760 VW dVAVy/' 

Furthermore, from the consistency conditions (see [ 6 ] for the meaning of the various terms) 

= ( 6 . 6 a) 

63 = (-F> + dih3 - 9/6i 4 + Tf - ( 6 . 6 b) 

it is clear that at two loops = 63 ^^ = 0. This, in conjunction with (13.4p . (16.21) and (16.41) . 
implies that = 0. These results have been verified by our explicit computations (I4.14p . Now, 
at two loops we can use (I4.14p and (j4.16p to obtain 


r(2) 


1 


1 


^ “ (647r3)2 1080 

and so using (16.61) we can compute 


91 ■ 


^( 2 ) _ l(2 ) _ g 

"13 ~ "14 ~ ’ 


^ 6 ( 2 ) ^ 


1 


1 


R3) ^ 1,(3) ^ __ 

^ 6 3 (647r3)3 6480 


(6.7) 


( 6 . 8 ) 
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With these results and using (13.41) with ()5.6I) . (I6.2p . and (16.411 we find that the three-loop contri¬ 
bution to o is 

= (64^3)3 64^0 (©" i@) • 

This shows that, just like d, a increases in the flow out of the trivial UV fixed point in our theory. 

There is one comment to be made about the value of the result in ()6.5p . It is a scheme- 
dependent quantity, in the sense that it is only defined modulo terms that are “exact” in the 
cohomology generated by the Weyl transformations — A^, i.e. up to local additions to the 
original action whose variations shift quantities in (13.6p . However, as shown in [U] in analogy 
with [5] these shifts are of the form 6d = zij/3 ^for zjj an arbitrary regular symmetric function 
of the couplings. Hence, at lowest order, and using (15.6p . we have dd ~ O(g^) which cannot 
possibly upset the conclusions of this section in perturbation theory. Moreover, equation (|3.6I) 
is of course unchanged by such shifts and in this sense is an invariant of the associated Weyl 
cohomology. 
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Appendix A. Conventions and basis tensors 

In this work we define the Riemann tensor via 

= (A-1) 

and the Ricci tensor and Ricci scalar as R^i, = and R = We also commonly use 

the Weyl tensor defined in d > 3 by 

kU/ii/po- R^ii/pa T (^_2 inp,[a^p]v T 'lv[p^a]p) T 1p[pda]u^ • 
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At mass dimension four we use the tensors 


^4 = id-2Kd-3) + R^) , 

F = , 

= (^-2K^-3) ^4^^^ _ 4(^ _ + 8//3M. + 8^4^. - , (^-3) 

H 211 U d—l RR/J-i' 1 RSfiu — Rff Rpu ) Hipv — R^ Rppcri/ ) 

H^pu = '^‘^Rpu , Hqpu = -^VpdyR. 

A complete basis of scalar dimension-six curvature terms consists of [36] 

Ki = R\ K 2 = RR^''Rpu, Ks = RR'^’^P^ Rp,p^ , K^ = R^^''R.pRP^, 

K 5 = RP'^RP^Rpp^u , = RP-'Rpp^rRf^^ , Kj = RP'^p^R p^r^R^^pu , 

Ks = RP''P^RrupcoR;‘^a, Kg = RV^R, Kw = R^^V^Rp., Kii = RP''P^V‘^Rp,p„, 

K ,2 = R^'^VpduR , Ki3 = VpR^p , i^i4 = VPRJ'P V,Rpp, 

Alls = VPR^P^^VpR.p^r , Alie = , K^j = {V^fR. 

In d = 6 a convenient basis is given by 

h = - ij ^2 + |ji^3 + - 1^5 - IK, + K, , 

h = 4^1 - %K 2 + ^i^3 + |i^4 - 1^5 - 3i^6 + Kj , 

/3 = + %K 2 - |d^3 - + 6 X 5 + 2i^7 - 8/^8 

+ f-K’g — SATio + 6 X 11 + 3Ari3 — GAT^ + SA'is , 

= K\ — 127^2 “ 1 “ 3-K^3 -|- 16-ftr4 — 24ii^5 — 24i^g -|- AK^ -\~ S-ftTg , 

Ji = GA'e — 3 Ar 7 + 12Ar8 -I- A'lo — TA'n — IIA '13 + 12 A'i 4 — dATis , (A.4) 

J 2 = -\Kg + Kig + lKi 2 + Ki3 , J^ = K^ + K^- ^Kg + |Ki 2 + ATu , 

Ji = —\Kg + ATii + |A'i 2 -|- Alls , J 5 = ATie , Je = Kn , 

Li = --^Ki + iAr 2 - Kq, Lg = -^Ki + ^Ar 2 , 

R 3 = — 0000 -^1 + 150 -^2 — 75 ATo + 4q ATs + Yg A'e , A 4 = — ATi + ^ A '3 , 

R 5 30 300 J wR'^ ’ -^15 ) 

where the first three transform covariantly under Weyl variations, and Eq is the Euler term in 
d = 6 . The J’s are trivial anomalies in a six-dimensional GET defined in curved space, and the 
first six L’s are constructed based on the relation 6 ^ J d®x ^/7 Ai^..._e = f d^x-yj^ 

In this paper we use the above basis for dimension-six curvature scalars, but, although it is 
not necessary, we define di, 2,3 in general d, because of our use of dimensional regularization. More 
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specifically, we define 


TUJ 

Tupuj *'*' fi a 


d‘^ + d-4 
~ (d-l)2(d-2)3 


3(d^ + d-4) 
(d-l)(d-2)3^2 + 


2(d-l)(d-2) 

3d 




{d - 2)3 


{d-2y 


:K^- 


d-2 


K& + K^, 


TUJ 

paruj *' •' pv 


h = W^'^P^Wp^ru^W 

8(2d-3) 

~ (d-l)2(d-2)3 


24(2d-3) 

^ (d-l)(d-2)3^2 + 




(d-l)(d-2) 

24 


(d-2)3 


(d- 2 y 


:K^- 


12 

d- 2 


Kq + Kj, 


16 


4d 




+ - ^^\WP'^P^Wp,p^) 

2{d^ + d + 2) 2(^2 + 13d-6) 2(d-3)^^ 4(d-10)^^ 

■-^1 + ~r, -TT71-- -, -^—4i3 + — —^^4 


(d-l) 2 (d- 2)2 


(d- l)(d- 2)2 


d- 1 


(d-2)3 


•“ - • K. , !B - M. - i'-?“.-» f. 


d- 2 


d- 2 


d- 2 


+ I (d " 2)(d - 3 )Ki5 + 2td-l)td-%^^^ ■ 


(A.5) 


These satisfy ( 5 crA, 2,3 = 6 (t/i_ 2,3 for any d for which they can be dehned. 


Appendix B. Coincident limits 

Here we collect the coincident limits a;' ^ a: of the Seeley-DeWitt coefficients an,ij{x,x') of (j2.22ji 

1 /2 

and the various functions (i.e. a{x,x') and Ayj^(x,x')) needed therein to solve the recursion 
relation (j2.25p . Most of these results can be found in [19ll2lj^ . and [T7j, though in these works 
only the single coupling case was considered. 

The fundamental quantities of interest on a curved background are the geodetic interval 
a{x,x'), whose “equation of motion” is m 

^dPadpa = a, (B.l) 


^^Reference m uses the opposite curvature convention as the one used in this work, (EH). This must be taken 
into account when comparing expressions, as odd powers of curvature will have an extra relative minus sign. 
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and the van Vleck-Morette determinant Avm( 3 ;,x') , which describes the rate at which geodesics 
coming from a point separate, follows from a corollary of the above equation for a{x,x'): 

VV + 2d^ad^AU^ = , (B.2) 

1 /2 

with d the spacetime dimension. Here we have abbreviated a{x,x') and Ayj^(x,x') as a and 
1/2 

Ayjy, respectively, and will continue to do so throughout the rest of the appendix. 

1 /2 

From these two equations we can construct the coincident limits of a and Ayjy. We will denote 
the coincident limit of a quantity X as [A] for brevity. Throughout the rest of this appendix, to 
make the coincident limits as concise as possible, we use the semicolon notation for the covariant 
derivatives, e.g. Now, since a measures a distance, we clearly have 

H=0. (B.3) 


Now using (IB.ip and differentiating as many times as needed, we obtain the following limits: 


O'. 


<7 




/ij — 0 , 1 

\^upa\ — T dipa-up) i 


kr- 


pup. 


= 0 , 




+ Ru 


= ^RuoR^u + ^RP^Rn 

4 / 


+ T^a + T-^p), 

4 p po-T p 
— TE-tin Uv 


_ 2 o P 

K J- ^n.ij \ n K . I 


W 7 pua^] = k;p^.V] + liR^'^Rppau " Rp^Rpu) , 


<7 


P 0-1 — 


;pup cr 


= \cr 


;p pi/a 


^ - liRP^Rppau - Rp^Rpu) - 


pu-,p 


P + 2R 


;pu ) 


= uiR^^Rpu - iR^'^p^Rpup. - 2R,^n , 

r p, V p (T'\ 152 176 i 296 44 tv^|80 7>"_i_3t^_i_ 148 

u p a \ — 315 -N4 315 -N5 + 315-^6 63-^7 + 53-^8 + 7-N9 + 105-^10 

- ^ All - 1^12 + f Ai3 + |Ai 4 - f Ai5 - ^Ai 6 - f Aiy . 

-1 /Q 

The coincident limits of derivatives on Ayyj follow from (|B.2p : 


[^Vm] “ ^ 


~ ^Rpu 1 


[A 


1/2 

VM;pi^p 


] — 12 {Rpu-,p + Rpp\u Rup;p) ) 




ppaiy 


+ Rp^^ Ry par) + 


3 R. 


+ 


20 


pu 


[^VM;pi'p^] ~ [^VM;p^pi'] + \ {RpRpu R^ Rppau) , 


1/2 


[A 


1/2 p u pi _ _\_ 


VM;p u p 


~ 216 *^1 60 + 60 R3 + 


60 


60 - 


2 K, + i,K,-RK, + ^Kr- 


189 


20 

189 


20 P’^’P ’ 


^8 -ids ^9 


— ^Rio + y All + yAi2 — ^ A13 — ^ A14 + ^ A15 + A Aig + -uRn ■ 


(B.4) 


(B.5) 


We have only listed the limits that are needed to compute the coincident limits of the Seeley- 
DeWitt coefficients a^^ij up to 03^^. We also note that we have explicitly checked all limits with 
those found in the aforementioned references and find agreement. 
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Another quantity of interest, which is indirectly related to our calculations in this paper 
through the details laid out in appendix [Cl is y = Knowledge of its coincident 

limits is necessary for the computation of (14.8|] . We find 




■301 


+ \R,^P , 

._ ]^( -ni^p p _ ryiypo-r p 'l _L 1 7? ^ 

15 ^^p]p ^ ^upaT]p) “T ^ , 

[y,tJ.'^iy ] = ^-^4 + — ^Kq + ^Kj — ^Ks — i^Aig ■ 

+ yA'ii + ^Ari2 — ^A'ls — + ^^17 ■ 


ToKio 


(B.6) 


Now we may proceed to the quantities that are directly related to the central computations of 
this paper, the Seeley-DeWitt coefficients an,ij{x,x') that characterize the propagator’s response 
to the curved background with metric 'jpuix)- Restating its fundamental and defining condition, 
(I2.25p . we have 

7l0,n,ij + d^(T dpQn^ij = ~^VM ^ik{^Yy[P^n—l,kj') j 


with initial conditions 


dpCjdPao^ij = 0 and [oo,ij] = 6ij . 


(B.8) 


The limits of these coefficients depend on the elliptic differential operator of the form My = 

. In the case of the Lagrangian (|4.2p . we have 

4‘=4‘h 

Xij — rUy + 9ijk4^k T CijR ■ (^'9) 


+ X, 


with X — 
y, with Ay — Q^j 


(Unless explicitly stated we will take cj) to represent the background field for the sake of 
compressed notation.) Note that Ay is symmetric. Ay = Xji. Then, 

= Q^^ij ~ Xij , 

[5^fli^y] — ijdpi^-^Rdij Xij^ , 

[V^9,ai,y] = X{\{RP^Rpp^^ + RpP^^Rupar) - R/Rpu + f + | V^5,R)5y - i V^^.Ay , 

[dpV^ai,ij\ = X{R''P^^V^R,p^r - R^^^XpRup - ^Rp d^R + 3dpX^R)dij - i^^V^Ay , 

[X^dpai,i,] = [dpX^ai,ij] + i- Ay), 

[(V ) oyy] = ^Ae + ^Ay — jjgAg + j^Ag — ^Aio + ^An 

“ ^-^12 — ^Ai 3 — 305-^14 + i|o-^15 — 2510-^16 

+ ^A^"V^a,Ay + XdPRdpXij - i(V2)2Ay . 

(B.IO) 
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for the relevant limits of For 02 ,we have 


[0‘2,ij, 

[d^CL2,ij_ 


[V'aa, 




= + 6V2ii),5,,- - i(i? + V^)Xij + ^XikXi, , 

= + iRd^R + ^df,V^R)6ij 

- ^d^iRXij) - ^d^X^Xij + li^Xikdf^Xkj + d^X^kXkj), 

= — 540-^3 — 1890-^4 “ ^-^5 + ^-^6 “ 3780-^7 + Tsg-^S “ T008-^9 ^TO-^IO 

“ Tio-^ll “ 140-^12 + + 420-^14 “ slo-^lS “ SMO^^ie “ 

- - R>^‘'R^,u + ev^i? + ddf^Rdp + 2,R^^'^Vpd^ + 5RV^ + |(v2)2)Xy 

+ lilXikV^Xkj + V^X.fc Xkj + d>^X,kdpXkj). 

(B.ll) 


Finally for 03 ^ jj we have 

M = ^{fK, - fK2 + + |i^4 + fiFs - f i^6 + fKj - fKs + llKg 

— 8 iFio + 12 iFii + 12 iFi 2 — 2 iFi 3 — 4iFi4 + 9Ki^ + ^iFie + ^^Kij^dij 

171/ TP TPP^ TP _L 5 /?2l 

30 Uv-^ il^ypo- R Rpu + 2 ^ ) (B 12) 

+ X^R + d^^Rdp + ^RX^ + iR^^'^Xpd^ + \{X^f)X,j 

+ ^{^RXikXkj + XikX'^Xkj + X'^XikXkj + d^XikdpXkj) 


^Xij^XkiXij . 

It should be noted that all derivatives are evaluated inside the coincident limits, i.e. the 
coincident limits of derivatives on the quantities are to be taken, not the derivatives on the 
coincident limits of said quantities. Actually, it is not difficult to convert between the two—the 
relevant equation was given by Christensen m and was, in fact, necessary for the computation 
of (I4.9|) . 


Appendix C. Coincident limits and divergences of prodncts of propagators 


In this appendix we give the e poles of the products of the propagator pieces G„(x,x') and 
Rn{x,x') found in (I2.3ip . As noted there, the full propagator in (12.150 is regular for any d at 
separate spacetime points. Explicitly, 


Go(x,x') 

1—1 

1 

"IS 

l-H 

1 

^ 1/2 

^VM 

^^d/2 

(2(t)'^/2-i ’ 

Gi{x,x') 

F(id-2) 

aV2 

^VM 

167r'^/2 


Rn{x,x') 

/ 1 


2 (-1)”-V"" 

e TT^{n — 2)! 



(C.l) 
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for n = 2,3. 

However, upon taking products with other Green’s functions, as in (14.81) . short-distance sin¬ 
gularities will arise that will have poles in e, in accordance with (IC.lh . The associated relations 
centrally depend on the coincident limit of inverse powers of a in non-integer dimensions. 


1 

(2cr(x, x')) 


^-5 2 vr '^/2 


5‘^{x,x '), 


(C.2) 


which is valid up to finite contributions as x' —>■ x. Here 5 oc e and the factor is inserted to 
preserve dimensions. This dependence can be seen from the a dependence in equations in (IC.ip . 
Varying powers of cr will arise depending on the product of propagators taken. A useful recursion 
relation can be obtained by differentiation and the use of (jB.ll) and (|B.2I) . It reads 


(V^ - Y) 



aP 


p{2p + 2- d) 


^1/2 
^VM 
crP+1 ■ 


(C.3) 


Equation (1C.311 can be used to obtain the poles in products of propagators. For example, if we 
multiply (1C.21) with Ayj^, act with — V and use (|C.3I) . we obtain 


^1/2 


Ai 


-5 


TT 


d/2 


(2^)i('=*-5)+i d dr(id) 


(V^ - Y)6° 


(C.4) 


which is necessary for determining the e poles in {GoGiGi){x,x') in d = 6 — e. 
Using these methods we can now list the relevant products as in m- 


(GoGo)(x,x') ~ , 

iGoG,)ix,x')^-^-J‘^{x,x'), 

(GoGiGi)(x,x') ~ (V^ + IR) 6\x,x') , 

iGoGiR 2 )ix,x') ~ 5‘^{x,x '), 

{GoGoRs){x,x') ~ (1 + n) 6\x,x'), 

iGoGoR 2 ){x,x') ~ I (V^ + IR) 5\x,x') 


(C.5) 


where ~ indicates that only the e poles are considered on the right side. In these expressions 
1 /2 

powers of Ayj^ that appear in the propagator products are commuted through to the delta 
function and then we use Ay^{x,x')6'^{x,x') = 5'^{x,x'). 


28 




























For our purposes we also need to know the divergent behavior of the products {GqGqGi){x,x') 
and {GoGoGo){x,x'). Their computation is performed after taking advantage of the fact that 
they appear under x, x^integrals, and thus we can integrate by parts at will. For example, for 
{GoGoGi){x,x') we need to hnd the poles in From (|2.3ip and (|C.4p . ()C.3p 

we see that we need 

Avm(V2 - Y f6‘^ = Avm + Y'^6'^] . (C.6) 

For the contribution of {GoGoGi){x,x') to (j4.8p we also have giki{x)gjki{x')ai^ij{x,x') under the 
X, x'-integrals. Following the procedure that led to (jC.Sh we would now commute Avm through 
the derivatives to the delta function. This is rather tedious, so we choose here to integrate the 
derivatives in ()C.6p by parts instead. This way we are be able to do the x^integral which will 
force the coincident limits of the various contributions that arise. The necessary results are then 
found in (IB.Sh . ()B.6I) and (IB.101) . 
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